Spin diode based on a single-wall carbon nanotube 
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Abstract 

Electronic transport through a single-wall metallic carbon nanotube weakly coupled to one fer- 
romagnetic and one nonmagnetic lead is analyzed in the sequential tunneling limit. It is shown 
that both the spin and charge currents flowing through such systems are highly asymmetric with 
respect to the bias reversal. As a consequence, nanotubes coupled to one nonmagnetic and one 
ferromagnetic lead can be effectively used as spin diodes whose functionality can be additionally 
controlled by a gate voltage. 

PACS numbers: 72.25.Mk, 73.63.Kv, 85.75.-d, 73.23.Hk 



Electronic address: weymann@amu.edu.pl 



Since their discovery, [jj carbon nanotubes (CNTs) have been extensively studied from 
both fundamental and application points of view, p, [J Owing to extreme flexibility of 
nanotubes (they can be either metallic or semiconducting), they have turned out to be ideal 
natural systems to study one-dimensional electronic transport in various transport regimes. 

Transport characteristics of a CNT contacted to metallic leads depend on the strength 
of CNT-leads coupling. Q| For very good coupling, transport reveals features typical for 
electron waveguides. In turn, for weak coupling between the nanotube and leads, CNT be- 

laves as a large quantum dot with many orbital levels participating in electronic transport. 

5I In the intermediate case and at sufficiently low temperatures, transport reveals features 
characteristic of the Kondo phenomenon. jf| When the leads are additionally ferromagnetic, 



transport properties of a CNT depend on the relative orientation of t 



re leads' magnetic mo- 



ments, leading to the so-called spin- valve effect, say , [nj, [ll|, Q, Q] Very recently 



Merchant and Markovic 



15| observed strong diode-like behavior in transport characteristics 



of a CNT coupled to two metallic leads, one being ferromagnetic (Co) and the second one 
nonmagnetic (Nb). Motivated by this experiment, in this Letter we analyze theoretically 
transport through a single-wall metallic carbon nanotube weakly coupled to one ferromag- 
netic and one nonmagnetic lead, and show that indeed such systems reveal features which 
are typical of spin diodes. By considering transport in the sequential tunneling regime, we 
show that both the spin and charge currents become strongly asymmetric with respect to 
the bias reversal. The magnitude and sign of this asymmetry depends on the bias voltage, 
and can be additionally controlled by a gate voltage, which is of particular interest from 
the application point of view. It is also worth noting that a similar spin diode behavior has 
been predicted for single semiconducting quantum dots. Q 

To describe the spin-diode features in transport characteristics, we consider the system 
consisting of a single-wall metallic CNT which is weakly coupled to one nonmagnetic (left) 
and one ferromagnetic (right) lead. Hamiltonian of the system has the general form, H = 
Hi + H-r + i^cNT + Hi. The first two terms describe noninteracting electrons in the leads, 
H r = £rko-<4 ko .c r k<T, for the left (r = L) and right (r = R) leads, with e r k a being the 
energy of an electron with the wave vector k and spin a in the lead r. To describe the CNT 
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in the weak coupling limit we employ the model Hamiltonian introduced in Ref. 18 



where N = E^a^wffi an d n nja = ^jcAni*- The discrete energy level in the subband ji 
(ji = 1,2), e^j, is given by = jA+(p—l)5, where A is the spacing between levels following 
from quantization in a particular subband, and 5 is the energy mismatch between the level 
sets corresponding to the two subbands. The charging energy of the nanotube is denoted 
by U and Nq is the charge induced by the gate voltage. The additional Coulomb energy of 
two electrons in the same level is described by 8U, while J is the exchange parameter. 

The tunneling Hamiltonian reads H T = Er=L,R Ek Y.^ai t rj c l\,a d M^+Kj d \ja c r^a), where 
t r j is the tunnel matrix element between the lead r and the jth level. Coupling of the 
jth level to external leads can be described by rjT- = 2n\t T j\ 2 p", where p a r is the density 
of states in the lead r for spin a. Defining the spin polarization p of the ferromagnetic 
electrode as p = (p^ — p^) / (p^ + p^), one can write the coupling of CNT to the right lead as 
r^j~' ) = a(l±p)T/2 for the majority (minority) electrons, while coupling to the nonmagnetic 
left lead as rjK = rj~j = T/2 (for all levels j). Since CNT may be coupled to electrodes 
with different coupling strengths, we introduce an asymmetry factor a, and assume a = 0.2 
(coupling of CNT to the ferromagnetic electrode is weaker than to the nonmagnetic one). 

When the CNT is weakly coupled to external leads, current flows due to the sequen- 
tial tunneling processes, except the Coulomb blockade regions, where cotunneling processes 
dominate over the sequential ones. However, both sequential and cotunneling contribu- 
tions to the current in the Coulomb blockade are small in comparison to the sequential 
current out of the blockade regime. The latter regime is of particular interest from the 
application point of view and is relevant to many recent experiments. Therefore, we re- 
strict the following considerations to the sequential transport regime, and use the real 
time diagrammatic method limited to the first-order expansion. Tunneling rates in the 
sequential approximation are given by the usual Fermi golden rule. In turn, the occupa- 
tion probabilities can be calculated from the master equation (WP) = T5 XX0 , [jj], 
where P is the vector containing the occupation probabilities, and \x) is a many-body 
state of the CNT. The elements of matrix W are given by W xx > = W^ x , + H 7 ^,,, with 



W XX' = 271 S CT Pr{fr(£ X - £ X') | E M j %j(x\dhjM\ 2 + I 1 ~ fr(£ X '- £ x)} | E A y trj(x\ d K>M\ 2 }, 

for X ^ X', and W r xx = -£ xYx W£ x , where / r (e) = i/^-MAbT + ^ and ^ ig the 
electrochemical potential in the lead r. The first (second) term in the bracket describes 
tunneling to (from) the nanotube from (to) the lead r. Apart from this, one arbitrary 
row Xo of matrix W has been replaced by (r, . . . , T) due to the normalization condition 
Tr{P} = 1. The sequential current flowin g th rough the nanotube can be then calculated 
from the formula I = e/(2n)Tr{W I P}, HQ where the elements of matrix W are given 
by W\ x , = [Q(N X , - N x ) - Q(N X - N x >)] (W* x , - W^), with N x being the number of elec- 
trons in state \x) and Q(x) denoting the Heaviside function. 

In order to model the single-wall metallic carbon nanotube we have taken the parameters 

n 

derived from the experiments of W. Liang et al. [5] We have also introduced a conversion 
factor, x, which relates the gate voltage to the electrochemical potential shift. 20] In Fig. [1] 
we show density plots of the absolute value of the current /, differential conductance G, 
absolute value of the spin current Is, and the differential spin conductance Gs, calculated as 
a function of the gate and bias voltages. The spin current is defined as Is = If — where Jo- 
is the current flowing in the spin-cr channel. Is is related to the angular momentum current 
Is by the relation I s = (h/2)I s /e. Due to periodicity, in Fig. [TJwe show only one sequence of 
the four-fold shell filling structure of single-wall CNTs. First of all, the Coulomb blockade 
regions are clearly visible - see the large and small black diamonds in Fig. [T](a). We also 
point on the inverse symmetry of the transport characteristics shown in Fig. [1] with respect 
to the center of the large blockade diamond. 

In Fig. [2] we show the bias voltage dependence of the transport characteristics (charge and 
spin currents and the corresponding differential conductances) for different values of the gate 
voltage, V g = —0.032 V and V g = —0.06 V, respectively. These characteristics correspond 
to the vertical cross-sections of Fig. [T] through the center of the first small diamond to the 
right of the large one, and to the resonance point between the large and small diamonds, 
respectively, see the dotted lines in Fig. [TJa). We recall that the large diamond corresponds 
to the situation when the sequence of the four electronic levels of the CNT becomes filled 
with four electrons, while the first small blockade diamond to the right of the large one 
corresponds to the situation when these four-level sequence is occupied by three electrons 
only. In turn, the first small diamond to the left of the large one describes the situation 
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when the next four-fold level-sequence is occupied by a single electron. 

As can be seen in the figures, the current flowing through the system is not symmetric 
with respect to the bias reversal. This is associated with the asymmetry in bare tunneling 
matrix elements for the spin- majority and spin- minority electrons between the CNT and 
ferromagnetic lead. Due to this spin-dependence, tunneling processes involving the spin- 
majority electrons are faster than those involving the spin-minority ones. On the other 
hand, the matrix elements for tunneling between the nonmagnetic lead and the nanotube 
do not depend on the spin orientation. Due to the spin dependence of tunneling processes, 
charge current is associated with a nonzero spin current Is = If — ij. flowing through the 
system. Spin current, similarly as the charge current, also reveals features typical of spin- 
diode, see Fig. [2](b) and (d). Thus, carbon nanotubes when coupled to one nonmagnetic and 
one ferromagnetic lead can be used as spin diodes, and their functionality can be controlled 
additionally by a gate voltage, see Fig. [TJ 

Let us now look more carefully at the two various situations shown in Fig. [2j In Fig.[2]^a,b) 
we show the bias dependence of the transport characteristics for the gate voltage correspond- 
ing to the Coulomb blockade region with three electrons in the last level-sequence of the 
CNT. We recall, that for the parameters assumed, electrons tunnel much easier through the 
barrier between the CNT and nonmagnetic lead than through the barrier between CNT and 
ferromagnetic lead. At small bias the system is in the blockade regime and current starts to 
flow (for both bias polarization) when the bias voltage exceeds a threshold. 20j For positive 
voltage electrons tunnel from the ferromagnetic lead to the nonmagnetic one. In Fig. [2] for 
V g = —0.032 V the step in current for positive bias is smaller than that for negative one. 
This can be explained by taking into account spin asymmetry in tunneling matrix elements, 
and difference in barriers between the left and right leads (barrier between the magnetic 
leads and CNT is larger). For positive bias transport goes mainly through the states | tli T) 
(lowest level doubly occupied and the next one occupied by a spin-up electron) and | fj,; 0) 
(the lowest level doubly occupied and the next one empty). In other words, an electron 
leaves first the CNT and tunnels to the nonmagnetic lead, and then a spin-up electron from 
the magnetic lead tunnels to the CNT. For negative bias, in turn, transport takes place via 
states | fl; Tl) (both levels of the sequence are doubly occupied) and | f J.; |) or | |; f J.) 
(one of the two levels is doubly occupied and the second one is occupied by a spin down 
electron). Now, an electron tunnels first to the CNT and then a spin up electron, either from 
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the lower or from the higher level, tunnels to the magnetic lead. Thus, there are now two 
channels for tunneling from the CNT to the magnetic lead, which makes the step in current 
at the threshold large for negative bias. For positive bias there was only one channel open 
for tunneling so the current step was smaller. When the bias increases further, additional 
channel becomes open also for positive bias and the currents become comparable. However, 
the two plateaus above the threshold for positive and negative bias are significantly different. 
In addition, the asymmetry in charge current leads to the corresponding asymmetry in spin 
currents, see Fig. [2£b). Moreover, the asymmetry for spin current is even more pronounced 
than for the charge current. As the charge current varies rather monotonically with the bias 
voltage (positive differential conductance), the spin current for positive bias drops with in- 
creasing voltage above the threshold (negative differential spin conductance). This behavior 
opens a functionality range for the system as a spin diode. 

Owing to the symmetry of the energy spectrum of the system, the situation is reversed 
when the gate voltage admits one electron on the last level-sequence (V g = —0.208 mV in 
Fig. [1]). Due to the corresponding particle-hole symmetry, transport for positive bias goes 
mainly through the states |0; 0) (empty levels) and | j; 0) or |0; f) (spin up electron on one 
of the two levels). For negative bias, on the other hand, transport takes place via states 
11; 0) (lower level doubly occupied) and | I; 0) (one spin-down electron on the lower level). 
Thus, now two channels open at the threshold voltage for tunneling from the magnetic lead 
to the CNT for positive bias, and therefore the corresponding step in current is larger than 
for negative bias. 

Consider now the situation shown in Fig. [2]^c,d). Since the system is now at resonance, 
there is no current blockade and the current increases immediately with applied voltage. 
Again both the charge and spin currents for positive bias are suppressed as compared to the 
currents flowing when the bias is negative. The arguments accounting for this bias-reversal 
asymmetry are similar as in the case of V g = —0.032 V. It is also worth to note, that in a 
broad voltage range the spin current for positive bias is much smaller than for negative bias, 
see Fig. E(d). 

We point that the spin diode behavior is observed for relatively low transport voltages. 
For higher voltages there are more states participating in transport and the current rectifi- 
cation is decreased. We also notice that the operation of the spin diode can be improved by 
increasing the spin polarization p of the right lead. For negative bias voltage when electrons 
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tunnel from nonmagnetic to ferromagnetic lead, once the nanotube becomes occupied by 
a spin-down electron the current is approximately given by I oc (1 — p)T/2, whereas for 
positive bias the current is / oc T/2. Thus, by increasing p, the rectification of the current 
is enhanced. An ideal CNT-based spin diode could be made by connecting the nanotube to 
half-metallic lead, where p — > 1 due to the energy gap in one of the two spin subbands. In 
addition, the operation of the diode can be tuned by sweeping the gate voltage, see Figs. [1] 
andd 

This work, as part of the European Science Foundation EUROCORES Programme SPIN- 
TRA, was supported by funds from the Ministry of Science and Higher Education as a 
research project in years 2006-2009 and the Foundation for Polish Science. 
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Figure captions 

Fig. 1. (Color online) The absolute values of the current I (a) and spin current Is = I^—Ii 
(c) in the units of Jo = eT/h, the differential conductance G (b) and the differential spin 
conductance G$ = dl^/dV (d), calculated as a function of the bias and gate voltages. The 
parameters are: A = 8.4 meV, U/A = 0.26, J / A = 0.12, 5U/A = 0.04, 5/A = 0.27, 
k B T/A = 0.025, p = 0.5, a = 0.2, x = 0.14, T = 0.2 meV, and I « 48.7 nA. 

Fig. 2. The current (solid line) and differential conductance (dashed line) (a,c), and the 
spin current (solid line) and spin differential conductance (dashed line) (b,d) as a function 
of the bias voltage for V g = —0.032 V (a,b) and V g = —0.06 V (c,d). The parameters 
are the same as in Fig. [TJ The insets illustrate the two subbands of the nanotube with 
the corresponding occupations - the CNT is occupied with iV = 3 electrons (a) [empty 
(filled) circles correspond to empty (occupied) charge states], while for V g = —0.06 V (c) 
the occupation number fluctuates between N = 3 and N = 4. 
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FIG. 1: (Color online) The absolute values of the current I (a) and spin current Is = I\ — ij, (c) 
in the units of Iq = eT/h, the differential conductance G (b) and the differential spin conductance 
Gs = dls/dV (d), calculated as a function of the bias and gate voltages. The parameters are: 
A = 8.4 meV, U/A = 0.26, J/A = 0.12, SUftM. = 0.04, S/A = 0.27, k B T/A = 0.025, p = 0.5, 
a = 0.2, x = 0.14, T = 0.2 meV, and I ta 48.7 nA. 
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FIG. 2: The current (solid line) and differential conductance (dashed line) (a,c), and the spin 
current (solid line) and spin differential conductance (dashed line) (b,d) as a function of the bias 
voltage for V g = —0.032 V (a,b) and V g = —0.06 V (c,d). The parameters are the same as in Fig. [TJ 
The insets illustrate the two subbands of the nanotube with the corresponding occupations - the 
CNT is occupied with N = 3 electrons (a) [empty (filled) circles correspond to empty (occupied) 
charge states], while for V g = —0.06 V (c) the occupation number fluctuates between N = 3 and 
N = 4. 
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